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As a model for nonideal behavior in the equation of state 
of QCD at high density, we consider cold quark matter in 
perturbation theory. To second order in the strong coupling 
constant, a„, the results depend sensitively on the choice of 
the renormalization mass scale. Certain choices of this scale 
correspond to a strongly first order chiral transition, and gen- 
erate quark stars with maximum masses and radii approxi- 
mately half that of ordinary neutron stars. At the center of 
these stars, quarks are essentially massless. 

PACS number(s): 12.38. Bx, 12.38. Mh, 26.60. +c, 97.60.Jd 



Strongly interacting matter under extreme conditions 
can reveal new phenomena in Quantum Chromodynam- 
ics (QCD). Compact stars serve as an excellent obser- 
vatory to probe QCD at large density, as their interior 
might be dense enough to allow for the presence of chi- 
rally symmetric quark matter, i.e., quark stars |l| pd| . 

The usual model used for quark stars is a bag model, 
with at most a correction ~ a s from perturbative QCD 
In the massless case, the first order correction cancels 
out in the equation of state, so that one ends up finally 
with a free gas of quarks modified only by a bag constant. 
If the bag constant is fit from hadronic phenomenology, 
then the gross features of quark stars are very similar to 
those expected for neutron stars: the maximum mass is 
« 2.M0, with a radius ps 10 km. 

In this article we consider quark stars, using the equa- 
tion of state for cold, dense QCD in perturbation theory 
to~ a 2 ||J|. These results are well known, and our only 
contribution is to use modern determinations of the run- 
ning of the QCD coupling constant ]l2| ]. At the outset, 
we stress that we do not suggest that the perturbative 
equation of state is a good approximation for the densi- 
ties of interest in quark stars. Rather, we use it merely 
as a model for the equation of state of QCD. 

To ~ a 2 s , there is significant sensitivity to the choice 
of the renormalization mass scale. Under our assump- 
tions, we find that this choice is tightly constrained by the 
physics. We consider two illustrative values of this pa- 
rameter. One choice corresponds to a weakly first order 
chiral transition (or no true phase transition), and gives 
maximum masses and radii very similar to that of neu- 
tron stars. The second choice corresponds to a strongly 
first order chiral transition |l3| , and generates two types 
of stars. One type has densities a few times that of nu- 
clear matter, and looks like the stars of a weakly first 



order chiral transition. In addition, however, there is a 
new class of star |7j,|i"o| , with densities much higher than 
that of nuclear matter . For this new class, the maximum 
mass is « I.Mq, with a radius ~ 5 km. Other models 
with nonideal behavior also generate small, dense quark 
stars [§]|. 

Assume that the chiral phase transition occurs at a 
chemical potential /i x [fl4"| . Our perturbative equation of 
state is applicable only in the chirally symmetric phase, 
when the quark chemical potential fj, > [i x . In this phase, 
the effects of a strange quark mass, m s w 100 MeV |fL6| , 
are small relative to the quark chemical potentials, fi > 
300 MeV. Thus we take three flavors of massless quarks 
with equal chemical potentials [l|Jl7||. 

The thermodynamic potential of a plasma of massless 
quarks and gluons was calculated perturbatively to ~ a 2 
by Freedman and McLerran Q and by Baluni g], using 
the momentum-space subtraction (MOM) scheme. The 
MOM coupling is related to that in the modified mimimal 
substraction scheme, MS, as J^Jlq |20| : 
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a s = g 2 /(47r), with g the QCD coupling constant, and 
A = 151/48- (5/18) Nf, with N f the number of massless 
flavors. (While we take Nf = 3, we give formulas for 



arbitary Nf.) In the MS scheme, to 
namic potential is then 
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where G = G - 0.536A/) + Nf In N f , G = 10.374 ± .13 
Iplfl , and A is the renormalization subtraction point. In 
MS scheme, the thermodynamic potential is manifestly 
gauge invariant. We take the scale dependence of the 
strong coupling constant, a s = a s (A) as |l2| , p2| : 
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u = ln(A 2 /Ay , ft = 11 - 
and 02 = 2857 - 5033 A^/9 



2JV//3, 0i = 51- 19JV//3, 
f 325^/27. The scale A^ 
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is fixed by requiring that a s = 0.3089 at A = 2 GeV JL2|; 
for N f =3,@ gives A^g = 365 MeV. 

All thermodynamic quantities follow consistently from 
Q(fx). The pressure is given by = — 0(/x), the quark 
number density by n(fj,) = (dp/d/i), and the energy den- 
sity by e = — p + iin. Given our stated assumptions, 
the only freedom we have in the model is the choice 
of the ratio A/fi. To illustrate this, we take the values 
A/ M = 1,2,3. 

For reasons which will become clear later, we find the 
choice A = 2fi especially interesting. Start with a very 
large chemical potential, such as fi — 100 GeV, for which 
a s ~ .095 (for the purposes of discussion, assume Nf = 3 
at this scale). The first order term decreases the ideal 
gas pressure by ~ 6%; the sum of the first and second 
order terms decrease the pressure by ~ 7% of the ideal 
gas value. Because the strong coupling constant runs 
relatively slowly with fj, at large /i, even at /i = 1 GeV, 
where a s ~ .31, the first order term decreases the ideal 
gas pressure only by ~ 20%; the first and second order 
terms, by ~ 30%. 

As can be seen from (|2|) , the perturbative expansion of 
the thermodynamic potential is an expansion in a power 
series not just in a s , but in a s log(a s ). The logarithm of 
a s arises from the plasmon effect, where the Debye mass 
squared m 2 D ~ a s /j, 2 . Because gluons at T = have four- 
dimensional phase space in loop integrals, however, the 
plasmon effect is relatively innocuous, and only produces 
logarithms, log(m£>//i) ~ log(a s ). 

This is in stark contrast to the perturbative expansion 
of the free energy at nonzero temperature, T^0. While 
there is again a plasmon effect, m 2 D ~ a s T 2 , because in 
loop integrals static gluons at T ^ have a three dimen- 
sional phase space, the perturbative expansion is not in 
a s , but in x/oj. The series in y/ct^ is much worse behaved 
than that at fi =/= 0, T = 0, and does not converge until 
very high temperatures [B3[ . The convergence appears to 
improve after resummation (24|,|n||2(J, or by using Pade 
methods ps| . 

Consequently, the perturbative series for the thermo- 
dynamic potential may be much better behaved at fj, =/= 
(and T = 0) than at T ^ |§6). This does not imply 
that a given value of a s , which is adequate to compute 
the thermodynamic potential, works equally well for all 
other quantities. In particular, the gaps for color su- 
perconductivity are nonperturbative, 4> ~ exp(— l/y/a^) 
p^ , [n| , and much smaller values of a s appear to be re- 
quired to reliably compute them p8|] . In QCD, effective 
models find that even when /i ~ 400 MeV, these gaps 
are at most ~ 100 MeV jl5|]. As the relative change in 
the thermodynamic potential is only ~ (</>//z) 2 , then, for 
the equation of state in QCD, color superconductivity is 
never a large effect. 

To truly know how well perturbation theory converges 
at fi 7^ 0, it is imperative to compute the thermodynamic 
potential to ~ a%. Unlike the case of T ^ 0, which is 
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FIG. 1. The total pressure, relative to the pressure of an 



ideal gas, po; including terms to order • 
as a function of fi; A = 2/i. 



a s and to order ~ a. 



sensitive to nonperturbative effects from static magnetic 
gluons from ~ on, at fi ^ (and T = 0), the entire 
power series in a s is well defined p7[ . 

We now use the perturbative calculation of the ther- 
modynamic potential for /i < 1 GeV. Since both terms 
~ a s and ~ a 2 have negative coefficients, as a s (ji) in- 
creases with decreasing [i, eventually the pressure van- 
ishes. While it is clearly invalid using perturbation the- 
ory when p = 0, it at least provides a well defined model 
of dense QCD. In Fig. [I] we show the pressure for A = 2/i; 
it vanishes at fj, c — 425 MeV, where a s ~ .65. This cor- 
responds to a quark density ~ 4.35po ; where po is the 
density of quarks in nuclear matter, ~3x .16/fm 3 . 

A weakness in our model is how to match the equa- 
tion of state for massless quarks, (|2|), onto that for mas- 
sive quarks and hadrons. The quark chemical potential 
must be larger than one third of the nucleon rest mass, 
minus one third the binding energy of nuclear matter, 
M > Mmin ~ 313 — 5 MeV. While the pressure vanishes 
at /x m in, hadronic (or quark) matter certainly exists, with 
nonzero pressure, for all /i > /x m i n . Thus we imagine that 
a very "soft" equation of state for massive quarks (and 
hadrons) matches onto the equation of state in (||) at 
some [i > [i c (see also below). Consequently, [i c cannot 
be much higher than /i mul . 

It is this which limits the choice of A/// in our model. 
For A = //, pb c = 767 MeV. It is absurd to think that the 
pressure of massive quarks could be small to densities 
~ 33po- Thus we do not consider this case further. 

For A = 3/i, n c = 300 MeV when a s - .6. By the 
Hugenholtz-van Hove theorem, when the pressure van- 
ishes, the ratio of the total energy to the baryon number 
i_s E/A = 3/x c . For iron, E/A = 930 MeV. Thus for 
A = 3/x, E/A = 3/x c = 900 MeV, and, as suggested by 
Bodmer and Witten strange quark matter is abso- 
lutely stable relative to hadronic matter. 

While possible, we prefer an alternate view. Our per- 
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turbative equation of state is valid only in the chirally 
symmetric phase, for fi > fi x . Perhaps when fj, < fi x , the 
true equation of state is close to our perturbative model, 
but vanishes smoothly as fi — > fi m - m . As discussed later, 
this is a model for a weakly first order (or no) chiral phase 
transition. 

The structure of a quark star is determined by the so- 
lution to the Tolman-Oppenheimer-Volkov (TOV) equa- 
tions [H]. For the TOV equations, all that matters is the 
relationship between pressure and energy density. This 
is shown in Fig. || for A//i — 2, 3. The numerical solution 
of the TOV equations gives the mass radius relationships 
of Fig. ||. For any solution to the TOV equation, the 
chemical potential reaches its maximum value at the cen- 
ter of the star; as one goes out in radius, the chemical 
potential decreases, and equals fi c (where the pressure 
vanishes) at the edge of the star. For A = 3/z, the maxi- 
mum mass is M max 2.14Mq. At this mass, the radius is 
Rmax ~ 12 km; the chemical potential at the center of the 
star is fi ~ 456 MeV, which corresponds to a quark den- 
sity of p max 4.87(00- When A = 2/i, M max w 1.05M©. 
At this mass, the radius is R max sw 5.81 km; the chemical 
potential at the center is fi » 649 MeV, corresponding to 
a quark density of p max ~ 14po- 

To help understand these results, it is useful to com- 
pare to the equation of state of a nonideal bag model: 

n(n)=-^a eSt i A + B eS ; (4) 

-B e ff is an effective bag constant, and the parameter a e e 
measures deviations from ideality. A common choice is to 
take a e ff from the thermodynamic potential to one loop 
order, with a fixed value of the coupling constant: a e s = 
1 — 2a s /ir ||. In a bag model, the relationship between 
pressure and energy density is linear, p = (e — 4_B)/3, 
irrespective of the value of a c ff . Thus we can uniformly 
scale p, e, and B together, so the maximum mass and 



radius satisfy a simple scaling relation, M max ~ i? max ~ 

i/s i/2 1- 

For the chemical potentials of relevance to a quark star, 
somewhat surprisingly we find numerically that the pres- 
sure in (0) can be well approximated by the effective bag 
model, (|I|). This can be seen from Fig. |2|, where the 
relationship between pressure and energy is very nearly 
linear. When A = 3/z, the pressure agrees with a bag 

1 /4 

model with B c ^ = 157 MeV and a c g = .626 to within 
2% for fi : 300 -> 470 MeV. This is close to the usual 
value in the MIT bag model, S 1 / 4 = 145 MeV @. 
When A = 2/i, the pressure agrees with a bag model 
with B^ 4 = 223 MeV and a cS = .628 to within 4% for 
fi : 425 -> 650 MeV. 

Consequently, the mass radius relationships for our 
quark stars agree well with a bag equation of state. To 
~ 5%, the maximum masses and radii scale according to 

1/2 

- l/BJx . The shape of the mass radius curve is also the 
same as for a bag model. Notably, light quark stars have 
small radii. This is because for light stars, M <C Mq, 
the chemical potential at the center of the star is near 
/i c , and the equation of state is controlled by that of 
massless fields, minus a bag constant. 

Our results for A = 2/z can be compared to other 
equations of state for dense QCD |8|-|l0||. All of these 
can be viewed as models in which there is nonidcality 
at a scale significantly higher than nuclear matter den- 
sities. Ref. H uses a Schwinger-Dyson model, and finds 
Mmax ~ O.7M and i? max w 9. km. Refs. §] and use 
models with massive quasiparticles, so that the masses 
act as a type of nonideality. Ref. || finds M max w 0.8M Q 
and i? max ~ 4. km; ref. Q finds Af max « 1.35 AT© and 
-Rmax ~ 10. km. We note that in Nambu-Jona-Lasino 
models, stars with a quark core do not arise, even at the 
maximum mass JTT{ . 

What about the manifestly nonperturbative phase in 
which chiral symmetry is spontaneously broken, /i < /i x ? 
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To understand this, consider an expansion in a large num- 
ber of colors The usual large N c limit is to let 
N c — > oo at fixed Nf. Since quark loops are suppressed 
in this limit, however, gluons are only affected by quarks 
when /i ~ N^ 4 or larger. This is in contrast to the 
transition at a nonzero temperature, which occurs at a 
temperature ~ N® We then consider a generalized 
large N c limit, in which N f -> oo at fixed N f /N c 
The quark thermodynamic potential is ~ — A^A^/i 4 , 
and the quark number density n ~ NfN c fi 3 . In this 
limit, (x x ~ Af?, as at nonzero temperature. 

For the purposes of our discussion, all that matters 
is that baryon masses wis ~ N c . The baryon chemical 
potential is related to the quark chemical potential as 
fiB = N c /i ~ AT C . When the baryons are nonrelativistic, 
so their Fermi momenta kf ~ 1, the baryon number den- 
sity is rig ~ cJ^fcy ■ The degeneracy of baryons is at least 
~ A 7 /, and could easily be larger, ~ N*. The baryon ther- 
modynamic potential is (naively) fis ~ — dsk^/mB- For 
fef ~ 1, however, everything is fine: the baryon density 
is ~ AT C that of quarks, and f2s < (l/A^fi. In terms of 
the quark chemical potential, however, (i = [j,b/N c ; with 
rrig = ms/Nc, fi w m g + fe? / '(2msA r c ) + . . .. That is, for 
fc/ ~ 1, the region in /i over which hadrons are a rea- 
sonable description is small, ~ 1/N%. Also, the binding 
energy of nuclear matter is automatically ~ l/N^. 

From this large N c argument, we conclude that a 
hadronic description is applicable only in a very narrow 
region of [i\ for larger fi, still of order one, a quark de- 
scription is appropriate. This need not be a true phase 
transition; rather, simply that the thermodynamic po- 
tential may be very difficult to compute in terms of 
hadrons, but relatively simple in terms of quarks. For 

1/2 

example, when kf ~ N c , so n — m q ~ 1/AT C , naively 

d B Nc /2 ~ N 1/2 n. This cannot be right — the 

thermodynamic potential of baryons cannot dominate 
that of quarks. The only resolution is that there are can- 
cellations — analogous to those which occur for baryon- 
meson couplings |35| ] — which greatly reduce the baryon 
thermodynamic potential, so that it is comparable to that 
of quarks. In other words, the hadronic thermodynamic 
potential must "soften" whenever fi — m q > 1/N%. 

Once one is away from this (narrow) hadronic window 
in fi, the appropriate equation of state for fi < \i x is that 
for massive quarks. There are then two possibilities. 

The first is that the thermodynamic potential for mass- 
less quarks matches, more or less smoothly, onto that of 
massive quarks. This requires either a weakly first order 
chiral transition, or perhaps just crossover. Below fi x , as 
f- Mmin the quark thermodynamic potential vanishes 
in a fashion typical of massive particles; within ~ I/A",? 
of Amin> a hadronic description is applicable. 

This is illustrated by the choice A = 3/i. For a star at 
its maximum mass, at the center /i s» 456 MeV; as the 
radius increases, \i decreases, until one enters a phase 



in which chiral symmetry is broken (assuming fj, x < 
456 MeV). Thus realistically, all of our "quark" stars 
have mantles with massive quarks, and then hadrons. For 
stars near the maximum mass, we assume this mantle is 
thin, and does not greatly alter its properties. As the 
mass decreases, however, the portion of the star in the 
chirally symmetric phase becomes small, until the entire 
star is composed entirely of massive quarks and hadrons. 
At this point, the relationship between the stars mass 
and radius is no longer like that of Fig. ||. Instead, it 
looks like that of nonrelativistic matter, for which the 
radius increases as the mass decreases. It is necessary for 
the chiral transition to be weakly first order, or a smooth 
crossover, for the mass-radius curve to be continuous. 

The second possibility is that the equation of state for 
massless quarks does not match smoothly onto that for 
massive quarks, with a strongly first order chiral transi- 
tion fisjj . As the thermodynamic potential approaches 
the ideal gas limit at large /i, and vanishes at /i m j n , 
this requires that the pressure is small at a value of 
fi x 3> fJ-min- Below fi x , by construction the pressure of 
massive quarks is small, with a "soft" equation of state 
& 

This occurs if A — 2fi. At the maximum mass, 
fi - 649 MeV at the center; if [i x < 649 MeV, most 
of the star is composed of massless quarks. As the mass 
of the star decreases, so does the amount in the chirally 
symmetric phase. If the chiral phase transition is strongly 
first order, eventually one jumps to a second branch, in 
which the chemical potentials are always < fi x through- 
out the star. Stars on this second branch are composed 
only of massive quarks and hadrons, with a maximum 
mass and radius like that of "ordinary" neutron stars. 
A strong first order chiral transition is necessary to en- 
sure that there are two, distinct branches. Using toy 
models for the thermodynamic potential, numerically we 
obtained solutions to the TOV equations which display 
two branches: we patched a thermodynamic potential for 
massive quarks, for fi < fi c w /i x , onto that for massless 
quarks, for \x > /i c . In this case, our stars of massless 
quarks constitute a third class of compact stars, after 
white dwarfs and "ordinary" neutron stars j?],[l0| . 

Most pulsars have masses - 1.5M Q [§. The MACHO 
project has also reported micro-lensing events for the 
Large Magellanic Cloud with masses M = 0.15-0.9Af Q 
p7f . For a weakly first order chiral phase transition, if 
MACHO events are hadronic stars, they must have large 
radii. For a strongly first order chiral phase transition 
]l3|] , MACHO events could be quark stars, with small 
radii, and pulsars might represent the second branch. We 
stress that our numbers for the maximum mass and ra- 
dius are meant only to be suggestive. Even so, we believe 
that our conclusions are qualitatively correct; a weak (or 
no) chiral phase transition leads to one type of compact 
objects, a strongly first order chiral transition, to two. 
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